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Abstract 



We consider the scalar delayed differential equation ex{t) = —x{t) + f{x{t — r)), where e > 0, r = r{x, e) 

l/^ • and / represents either a positive feedback df/dx > or a negative feedback df/dx < 0. When the delay 

P^ I is a constant, i.e. r{x, e) — 1, this equation admits metastable rapidly oscillating solutions that are 

pH ' transients whose duration is of order exp(c/e), for some c > 0. In this paper we investigate whether this 

metastable behavior persists when the delay r{x, e) depends non trivially on the state variable x. Our 

conclusion is that for negative feedback, the persistence of the metastable behavior depends only on the 

way r(x, e) depends on e and not on the feedback /. In contrast, for positive feedback, for metastable 

solutions to exist it is further required that the feedback / is an odd function and the delay r{x, e) is an 

^ I even function. Our analysis hinges upon the introduction of state dependent transtion layer equations 

l/^ ■ that describe the profiles of the transient oscillations. One novel result is that state dependent delays may 

lead to metastable dynamics in equations that cannot support such regimes when the delay is constant. 

^ ' Key words: metastability, state-dependent delayed differential equation, singular perturbation. 
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1 Introduction 

State dependent delay differential equations (DDEs) of the form 

dx 
^--j-{t) = ex{t) = —x{t) + f{x{t — r)), where r = r{x, e) = tq + f{x, e) , (1) 

appear as models in economics, physics and biology [H [251 [23l [Ml El 13] . Many studies have dealt with the 
asymptotic dynamics of these equations, establishing existence, stability and profiles of the so called slowly 
periodic oscillations when / represents a negative feedback (see, for instance, [U [31 [20l [22l |33l |43l [26l [311 ESI 
[38j). There are also results on existence and stability of periodic solutions, and on the convergence of most 
solutions to equilibria in the case of monotone positive feedback |211 141] . For a review on DDEs with state 
dependent delays see [19]. In contrast with these previous studies, in addition to the long term behavior of 
solutions, we are also interested in the transient dynamics of Eq. ([1]): we explore conditions under which the 
system displays metastable transient oscillations prior to convergence to its asymptotic state. Typically, in 
such a regime, the system engages in seemingly stable and sustained oscillations for a time of order exp(c/e) 
for some constant c, before eventually converging to a stable equilibrium or a stable periodic solution. 

It is well known that delayed feedback can lead to oscillations in a system that would otherwise remain 
at rest. Much focus has been on delay induced sustained undamped oscillations, such as periodic ones. 
However, the delay can also produce transient oscillations, that end up vanishing with time. While these 
had been reported in several studies (see for instance [42^ [TSl \T6\ [39] and references therein), for a long 
time they had not been considered as a topic of investigation on their own. In [39], the term DITO was 
coined for this phenomenon both as an acronym for "delay induced transient oscillation" and as a means to 
emphasize its repetitive nature. In the same paper, it was argued that long lasting DITOs could strongly 
alter the dynamics of the system and, for instance, interfere with information retrieval in neural networks. 
This work together with the ones cited above initiated intrinsic interest in the phenomenon. Since these 
early studies, DITOs have been reported in other systems, such as for instance |34[ l35j. In support of claims 
in [l2] and numerical explorations in [39], through analytical treatment of a specific system of DDEs with 
piecewise constant feedback, it was shown, in [30], that some DITOs could persist for durations of the order 
exp(c/e). Such long transients would outlast any observation window when e is small enough and would 
be undistinguishable from (nearly) periodic solutions. In analogy with long lasting transient oscillations in 
partial differential equations, DITOs with such exponential lifetimes are referred to as metastable [17]. 

Metastable oscillations have been thoroughly studied in many systems such as scalar partial differential 
equations [9l[l4]. However they have received far less attention in DDEs. Here we examine their occurrence 
in DDE ([1]). Our previous works characterized the conditions for metastability in equations with constant 
delays for both positive and negative feedbacks |17[ I18j . It further revealed a remarkable difference in 
metastability between the positive and negative feedback equations, that had escaped earlier reports and 
explorations. Indeed, a positive feedback DDE presents metastability only when the feedback / satisfies a 
special symmetry property that holds, for instance, when / is an odd function. This condition is absent in 
DDEs with negative feedback (see also [361 [SZ] for related results) . One consequence of this result is that 
not all DITOs are metastable. 

All these previous works deal with constant delays. As far as we know, the present paper is the first one 
to investigate and show the existence of metastability for state-dependent delay equations. Our goal in this 
paper is to study the transient dynamics of Eq. ([1]) when the parameter e > is small or, equivalently, rg 
is large. In previous singular perturbation analyzes of DDEs with state dependent delay [28l [29l [30], the 
small parameter e appears only in the left hand side of the equation. Here, the novelty resides in the fact 
that e appears in Eq. ([T|) in two different places: multiplying the time derivative of x, which characterizes a 
singular perturbation problem, and as an argument of f. One major contribution of this work is to highlight 
the fact that the dynamics of Eq. ([1]) strongly depend on the way f depends on e. 

This paper is organized as follows. In section 2 we present and discuss some previous results on the 



existence of metastable solutions in the context of DDEs with constant delays. These previous results indicate 
that there are two essential ingredients for the occurrence of nietastability. The first is the existence of rapidly 
oscillating periodic solutions (in a sense to be defined in section [2|) as the singular parameter e tends to 
zero. The second is that the shapes of these oscillations tend to stereotyped "square-wave" profiles as e ^ 0. 
These profiles are heteroclinic solutions of certain transition layer equations. A key consequence of this is 
that nietastability can be characterized through the analysis of these transition layer equations. In sections 
[3] and [4] we examine the same two ingredients namely the existence of rapidly oscillating periodic solutions 
and transition layer equations for state dependent delay equations. More precisely, in section [3l using a Hopf 
bifurcation theorem, we examine the existence of rapidly oscillating periodic solutions for Eq. ([T]), regardless 
of the choice of R and its scaling with the small parameter e (details of the computations are presented in the 
appendix A). Then we link the amplitude of rapid periodic oscillations to the possible existence or absence 
of nietastability in some of the state-dependent DDEs under study. In section H] we introduce transition 
layer equations to describe nietastability (more details are given in Appendix B), similarly to what has been 
done for DDEs with constant delay [101 [271 IIH] • We show that symmetry conditions on solutions of these 
equations can be used to characterize the condition on the feedback function /, and the delay function r, 
required for the existence of metastable oscillating solutions. One novel observation, specific to equations 
with state dependent delay, is that the existence of metastable solutions may depend on the behavior of the 
function r when e — )• 0. In section [5l we present systematic numerical explorations of the transient dynamics 
of Eq. ([1]). Section [6] contains a discussion and a summary of the results. 

Prior to section [21 we go over some assumptions and notations. In most applications the functions / and 
r are differentiable in x, here this is always assumed. Furthermore, in accordance with our previous results 
pTl [18] , we assume that / is either a monotonic increasing (positive feedback) or a monotonic decreasing 
(negative feedback) function satisfying the following hypotheses. 

Positive Feedback 

f'{x) > 0, /(O) = 0,/'(0) > 1, and there exist a > 0, 6 > 0, such that /(-a) = -a, f{h) = b, 
< f'{-a) < 1, < f'{b) < 1, and /(x) / x for x € (-a,0) U (0,6). 

Negative Feedback 

f'{x) < 0, /(O) = 0,/'(0) < -1, and there exist a > 0, 6 > 0, such that /(-o) = b, f{b) = -a, 
0< f'{-a)f'{b) < 1; |/(/(x))| >xforxG (-a,0)u(0, 6), and |/(/(x))| < x for x G (-oo, -a)U(6, oo). 

In the following, a positive or negative feedback / will be called symmetric if it is an odd function. Under 
these hypotheses, both the dynamics of the map / : R — t- M, and the DDE Eq. ([T|), with constant delay 
are well understood. In particular, the map / in positive feedback case has only three fixed points: x = 0, 
unstable, and x = —a, x = b, stable. For negative feedback, the map has a single unstable fixed point x = 0, 
and a single attracting 2-periodic orbit {—a,b}, R{x) = corresponding to constant delay. Without loss of 



assume that the delay function is 



generality, one can use the time scaling s = —, and consider Eq. ([TJ with rg = 1. In the following we shall 



r(x,e) = l + r/(e)i?(x), (2) 

where 77 is a non negative and smooth function on e > 0. In this way, different scalings between the delay r 
and the parameter e can be analyzed by changing the function rj. Typically R will be chosen as a smooth 
nonnegative function. 

2 Metastability in DDEs with constant delay 

Metastability is a concept that appears in several branches of physics and mathematics. Schematically, a 
state is called metastable if it is transient, that is, it is eventually transformed into another one, but this 
transformation is on such a slow time scale that it is not perceived in normal observation windows. One of 



the first systems of relevance to the present work, in which metastability had a full mathematical treatment, 
is that of scalar parabolic equations [9ill4j. Metastable solutions of DDEs ([TJ with constant delay r(x,e) = 1 
share a number of features with those of the partial differential equations |17t [T8] . In this section, we provide 
an overview of metastable solutions in scalar DDEs with constant delay, that will serve as a basis for the 
comparison and analysis of the case of state dependent delays. 

Throughout the remainder of this section, we refer only to DDEs ([1]) with constant delay set to one. 
Roughly speaking, metastable solutions of these equations are trajectories that evolve close to unstable 
manifolds of unstable periodic orbits, and they seem to be periodic, while in fact they are endowed with a 
drift of order exp(c/e). In section \2A\ we describe the relationship between rapidly oscillating periodic orbits 
and metastable oscillations, and in section [2^2] we show how transition layer equations allow to quantify the 
slow drift of such oscillations. 

2.1 Metastability and rapidly oscillating periodic orbits 

To make the description more concrete, consider, for instance, the positive feedback case f{x) = ^ arctan(5x). 
For any value of e > 0, Eq ([T]) has exactly three equilibria, x = —a, x = 0, and x = a, where a is the positive 
solution of X = iarctan(5x). For all values of e > 0, the equilibria x = —a and x = +a are locally 
asymptotically stable and x = is unstable. For all e sufficiently large the unstable manifold of x = is one- 
dimensional and as e decreases, x = undergoes an infinite number of Hopf bifurcations at ei > e2 > . . . > 
such that the dimension of the unstable manifold of x = at e G (e„, e„+i) becomes 2n + 1. Let x„ denote 
the periodic solution that bifurcates from x = at e = en- If Xn(i) = then Xn{t + s) has exactly 2n zeroes 
for s G [—1,0]. The branches of periodic orbits that appear at these successive Hopf bifurcations can be 
extended up to e = 0, the period converges to r = ro = 1 and the amplitude converges to a > [2\, and the 
oscillations tend to a square-wave-like shape when e — ?> 0. So, Eq. ([1]) admits periodic solutions with a large 
number of zeroes in a time interval of length one, provided that e > is sufficiently small. Given e > the 
global attractor of Eq. ([T|) consists of its set of equilibria and periodic orbits, and their finite dimensional 
unstable manifolds. These solutions are organized in a peculiar way: the global attractor admits a Morse 
decomposition, with each Morse set containing a periodic orbit or an equilibrium. The direction of the flow 
on this decomposition is such that the number of zeroes of solutions is a non increasing function of time 
[USl]. A similar description of the periodic orbits, their branches, their shape, and the organization of the 
trajectories on the global attractor holds for Eq. ([1]) with negative feedback |261l27j. 

Now that we have depicted the long term dynamics of DDEs with positive and negative feedback, we can 
describe the way metastable solutions appear in Eq. ([1]). For e > small, an initial condition ip : [—1,0] — )■ M 
to Eq. ([T]) with 2n zeroes gives rise to a solution Xi(-,(/7,e) : [—1,0] — >• M that after a time of order one is 
pointwise close to a function in the unstable manifold of x„, that has a square- wave-like shape close to that 
of Xn- The dynamics of x„ is approximately periodically oscillatory, but the slow motion of the solution 
along the unstable manifold of x„ eventually annihilates a pair of zeroes of the solution, which takes a time 
of order exp(c/e). After the annihilation of the two zeroes the solution drifts along the unstable manifold of 
another periodic solution Xn-i, and this process repeats itself until the solution eventually approaches one 
of the two stable equilibria x = —a or x = a for positive feedback, or a slowly oscillating periodic orbit for 
negative feedback. Metastability means here that the solutions have a (fast) oscillatory transient time that 
grows as exp(c/e) when e — )• 0. 

2.2 Metastability and transition layers 

Metastability not only depends on the existence of periodic solutions (i.e. the qualitative geometry of 
the phase portrait), but also on the quantitative dynamics along their unstable manifolds. Metastability 
happens only when the rapidly oscillating solutions are square- wave-like, in which case their jumps have been 
understood as transition layer phenomena, and described and analyzed using transition layer equations. 



Suppose that / is a positive feedback and that x{t) is an oscillatory metastable solution of ([T]), with 
r = 1, that jumps from x{t) sa 6 > 0, for t < 0, to x{t) ^ —a < 0, for t > 0, with x(0) = 0. The 
approximately periodic behavior of the metastable solution implies that x{t) w x{t + 1 + pe). When e is 
small, this and Eq. ([1]) imply 

flnr 

e^{t) = -x{t) + f{x{t - 1)} p. -x{t) + f{x{t + pe)}. 

Rescaling x{t) as (j)~{t) = x{et), the above equation for x{t) implies that <p~ must satisfy the transition layer 
equation 

r(t) = -r(i) + /(r(i + P")), (3) 

where p~ > is an unknown constant and the function (j)~ must satisfy the boundary conditions limt^_oo 4'~ (t) 
b and limj_>.oo 4'~i't) = —cl- The same ideas apply to a jump from x{t) Ri 6 > for t < 0, to x{t) f« —a < 
for t > 0, and leads to the existence of an increasing transition layer solution, i.e. a solution to 

(^+(t) = -(/>+ (t) + /(<A+(t + p+)) (4) 

where p~^ > is an unknown constant and the function 0^ must satisfy the boundary conditions limt->-co 0^(i) 
—a and limi_>.oo 4'~^{t) = b. 

The constants p"^ and p^ are drift velocities of ascending and descending sign-changes (zeroes) of an 
oscillatory solution. If an oscillatory solution x{t) of Eq. ([1]) satisfies 3;(t^) = and x'(t^) > 0, then one 
expects to find a time t^ ss t+ + 1 + ep^ such that x{t^) = and x'{t^) > 0; if x(t^) = and x'(t(^) < 0, 
then one expects to find a time t^ si t~ + 1 + ep~ such that x{t^) = and x'{t^) < 0. This suggest that 
the symmetry condition p"*" = p~ should be associated with metastability of oscillating solutions. Indeed, 
exponential duration of oscillatory transients is proven in [18] for scalar DDEs, with monotone positive 
feedback, by an estimate of the type \t^ — {to + 1 + ep^)| < e i so that when p+ = p^ oscillatory solutions 
are close to a 1 + e/9 periodic solution up to an exponential order. 

Finally, we discuss transition layer for negative feedback function. Transient oscillations are square- wave- 
like and they have an approximate period T ~ 2 + eC. If an oscillatory solution x(t) of Eq. ([T]) satisfies 
x{to) = and x'(to) > 0, then one expects to find a time ti ^ to + l + ep~^ such that x{ti) = and x'{ti) < 0, 
and then a time ^2 ~ ^i + 1 + ^P~ such that x{t2) = and x'{t2) > 0. So the transition layer equations for 
the increasing and decreasing transition layer solutions are coupled: 

P+(t) = -</,+ (f) + /(</>-[t + p-]), 

r{t) = -r(i) + /(<A+[i + P+]), ^ ^ 

where (f)'^ is increasing and (j)~ is decreasing on M, with limt_^„oo (p'^it) = ~Oi hmt-i-oo (p'^it) = b, limt_^_oo (p^it) 
b and lim^^^oo (^) = ~Oi with 0^(0) = 0^(0) = and p^ are unknown real constants. For negative feed- 
backs, given an oscillatory solution x{t), an (ascending) zero x(to) = with x'(to) > gives rise to a 
(descending) zero x{ti) = with x'{ti) > 0, and then to another ascending zero x(t2) = with x'{t2) > 0, 
with ti « to + 1 + ^P^ and ^2 ~ ^i + 1 + ep • So, in first order, the drift speeds of ascending and descending 
zeros are identical and equal to 2 + e(p^ + p^). Hence the symmetry condition, equivalent to the condition 
p^ = p~ in the positive feedback case, is here p'^ -\-p~ = p~ +p~^ , and it is obviously always true, irrespective 
oi p+ andp- [B]- 

In summary, the analysis using transient layer equations reveals and explains that metastability manifests 
itself in different ways in positive and negative feedback systems. Indeed, positive feedback DDEs present 
metastability only when a very special symmetry property of the transition layer problems ^ and dH is 
satisfied, which holds for instance when / is an odd function, while there is no such restriction in the case 
of negative feedback equations ([5]) [171 H] • 



3 Existence and amplitude of periodic solutions 

Based on literature results and novel results that follow in this section, we conjecture that, provided the delay 
satisfies a number of classical technical assumptions [281 EHl EQl [201 [HI [El [SI [7], the geometric organization 
of the phase portrait of DDEs with state dependent delays and monotone feedback is similar to that of 
equations with constant delays. Namely our three conjectures are: (i) Equation ([T|), with f(x,0) = 0, 
sustains branches of periodic solutions, that appear at the equilibrium by successive Hopf bifurcations 
and exist until e ^ 0, with amplitudes and periods that converge to some non zero limits, (ii) A Poincare- 
Bendixson like theorem holds, and as a consequence the global attractor of ^ is composed of equilibria, 
periodic solutions, and their unstable manifolds, (iii) The global attractor of ^ has a Morse decomposition, 
it is ordered by a discrete Lyapunov functional, and it is composed only of equilibria of ([1]), the periodic 
solutions described in the first conjecture, and connections from more-rapidly oscillating periodic solutions 
to less-rapidly oscillating periodic solutions. 

With few hypotheses on the feedback /, we have shown that a local Hopf-bifurcation theorem of [12] 
applies to ([T]). This gives an essential element in the proof of conjecture (i): for some decreasing sequence 
efc — >■ 0, a Hopf bifurcation occurs at the zero equilibrium of ([TJ each time that e crosses one e^, which gives 
rise to an oscillating periodic solution with 2k zeroes per period. The corresponding theorem is rigorously 
stated in section \3A\ and proved in appendix [Al Admitting the existence of branches corresponding to these 
periodic solutions, in section [3121 we show that if we have ry(e) = ce + o(e) in ([2]) for some < c < +oo, then 
the amplitudes of all periodic solutions converge to the same non zero limit when e — )• 0. On the contrary, 
following a proposition of [29], we show in section [331 that when r/(0) ^ 0, the amplitude of periodic solutions 
with 2k zeroes per period (along the branch that appears at e^), is bounded from above by some constant 
Ck (independent of e) such that Ck — > 0. 

fc— ^+oo 

3.1 Sequence of Hopf-Bifurcations for Eq. ([TD 

The PhD thesis of M. Eichmann [,,12j contains a local Hopf-bifurcation theorem for state-dependent DDEs 
which implies the following. 

Theorem 1. Suppose that f is C'^{'R,R), r :]0, l[xC°([-M,0],R) -^MisC^, andr :]0, l[xC^([-M,0],M) -^ 
R is C^, for M > 0. Suppose that /(O) = and |/'(0)| > 1. Then there is a decreasing sequence 
{€k)keN converging to zero, such that for any k > there is an open interval ] — r/fe,r/fe[ and C^ map- 
pings y* :] - Vk,Vk[-^ C^{[-M,0],M.), e* :] - r?fc,%HO, 1[ and w* :] - Vk,Vk[-^ ^, with y*{0) = x* = 0, 
e*(0) = €k and w*{0) = Pk = I™(-^A:) such that for any u €] — rj^^rikl there is a periodic solution to 

e*{u)x'{t) = -x{t) + f{x{t - r{e*{u), xt))) 

with initial condition xq = y*{u) and with frequency ^2n ■ 

We remark that the delay function r{e,xt) in the theorem above is more general than that in Eq. ([2]), 
and that the function / does not have to be positive or negative feedback. 

The theorem is proved in appendix [3 Here we first justify the existence of the critical values efc, and 
show that the associated frequencies y^ = f^ — > -|-oo. Suppose that /(O) =0 and that |/'(0)| > 1, and 

let X* = be the unstable steady solution of Eq. ([T]). The linearization of Eq. ([T]) at x* = is 

ey'it) = -y{t) + /'(O)y(t - ro) 

where ro = r(0) = r{x*) = 1. The characteristic equation associated to this linearized equation is 1 + eA = 
f'{Q)e~^ , or, equivalently, with A = a -|- i/3 G C 

1 + ea = /'(O)e-" cos(/3) 
ej3 = -/'(0)e-°sin^, 



This is the same characteristic equation as for the constant-delay equation, and there exists a sequence 
Efc — > such that for each e = e^ the characteristic equation has a single pair of solutions on the 

k—^ + OD 

imaginary axis A = ±i/3fc, /3fc > 0. Moreover, /3fc — )■ cxd as A; — >• oo. Since (3k is the angular frequency of the 
periodic orbit unfolded at e^, theorem [1] implies the existence of rapidly oscillating periodic solutions of Eq. 
([T]) as e tends to zero. 

3.2 Case 77(0) = with < r]'{0) < +00 

With the previous assumptions on t], one has r{x, e) = 1 + r]{e)R{x) ~ 1 + er]'{0)R{x) when e ^ 0. Without 
loss of generality, we assume ry'(O) = 1, so that the delay is r(x,e) = 1 + eR{x) in Eq. ([1]). 

Suppose that there is an eo for which Eq. ([T]) has a periodic solution xo{t) with period Tq. Given an 
integer n > 0, Eq. ([TJ and the periodicity of xo{t) imply that 

(^oxo{t) = -xo{t) + f{xo{t - 1 - nTo - eoR[xo{t)]}) . 

Then, rescaling time as t = t/(l + nTo), we obtain that x„(t) = xo[(l + nTo)f| satisfies the equation 

'° '^''" (t) = -Xnii) + / f X„ It - 1 - -J^R[xr,{i)] 



l + nTo dt ' ' ^ ' VI l + nTo 

Therefore, for e = e^ = eo/(l + "■^o) Eq. ([1]) admits the periodic solution x„(t) = xo[t(l + nT^)] with period 
Tn = Tq/{1 + nTo). It shows that the branches of rapidly oscillating periodic solutions can be obtained from 
the first branch of periodic solutions. Hence, assuming that the first branch exist up to e — >■ 0, it follows 
that all the other branches exist, and the amplitude of periodic solutions along all theses branches converge 
to the same positive limit when e — )• 0. 

In this sense, Eq. ([T]) sustains "large amplitude" rapidly oscillating periodic solutions as e tends to zero 
when the delay function is of the form r(e, x) = l + eR{x) or r(e, x) = l+er/'(0)-R(x). As mentioned in section 
[21 this is one of the signatures of the existence of metastable solutions in the case of DDEs with constant 
delay. So this strengthens the similarity of DDEs with state dependent delay and DDEs with constant delay, 
thus giving support to the possibility of metastability in the case r/(0) = with < rj'[{)) < +00. 

3.3 Case 77(0) 7^ 

The situation for r/(0) 7^ is different from the one depicted above, and this can be understood thanks to 
proposition 3.4 in [29], which precludes the existence of large amplitude rapidly oscillating periodic solutions 
as e — )• 0. First, we recall that proposition in [29], and then we discuss its consequences in terms of the 
amplitude of the periodic solutions. 

Proposition 1 (Mallet-Paret, Nussbaum). Suppose that the feedback f is C''(M,]R) and that the delay 
function r{e,x) is Lipschitz regular in x. Let x{t) satisfy equation ([1]) for t € M for some value of e, and 
suppose there exist oq < oi and to < ti < t2 < ts such that x{ti) < oq for even i, and x{ti) > oi for odd i, 
for < i < 3. Then 

max r(-,e) — min r(-,e) < 3(^3 — to). 

[a(),ai] [ao,ai] 

Proposition 2. Suppose that the feedback f is C'^(M,M) and that the delay function r(e,x) is Lipschitz 
regular in x. Let r{x,e) = 1 + rj{e)R{x), and suppose that R{x) is not constant on any interval (e.g. 
R{x) = r^x , k>l and r^ ^ Oj and r]{e) ~ 77(0) 7^ 0. 

Then, there is a function if, depending only on r]{0) and R, with (p{T) — > 0, such that for any x{t) be 

a periodic solution of Eq. ([7P with period T, and oq = minx(t) and oi = maxx(t), we have 

\ai -aol < i/3(r). 



In particular, if f and r satisfy additionally the hypotheses of theorem [71 the periodic solutions x^ that 
appear when e = e^ (see theoremU^ have periods Tj. — > and amplitudes max{xfc(t)} — min{xfc(t)} — > 0. 

fc— >oo fc— >oo 

Proof. To see why proposition [1] precludes the existence of large-amphtude rapidly oscillating periodic so- 
lutions, let x{t) be a periodic solution of Eq. ([T]) with period T, and let oq = minx(t) and oi = maxx(t). 
Choose to such that x(to) = oq, choose ti such that t^ < ti < t^ + T and x(ti) = oi, and t2 = to + T 
and ^3 = ^1 + T. As in section [H let r(x,e) = 1 + r]{e)R{x), and suppose that i?(x) is not constant on any 
interval, e.g. R{x) = r^x'', A; > 1, r^ / 0, and T/(e) '^ 7/(0) / 0. Then 6T > 3{ts — to) > r/(0)|rfc|Aij/3, where 
Ar = max^g[„g^„^] R{x) - min^g[„(,^„^] R{x) > 0. Given that 77(0) 7^ 0, we have |rfc|Aij/3 < ^. 

This estimate relates the period of oscillations indirectly to their amplitude (through A/j) : the faster 
the periodic oscillations are, i.e. the smaller T is, the smaller their amplitude is. The Hopf bifurcation at 
e = efc gives rise to a periodic solution with period T^ = 4- — > (see theorem [1] and appendix ffl, and the 

estimate above implies that so does their amplitudes oi — uq — > 0. D 

We argue that this result indicates that DDE ([T|) with state dependent delay, and 77(0) / 0, cannot 
support metastable transient oscillations that resemble those of DDEs with constant delays. This is further 
supported by the fact that the profiles of oscillations in this case are not solutions to usual transition 
layer equations (see discussion in next section U]) and confirmed through extensive numerical investigations 
(section [5]) . 

4 Transition layer and metastability 

In this section, we refine our previous analysis of the conditions under which DDE ([T|), with state dependent 
delay, can support metastable oscillations through the introduction of transition layer equations. Such 
equations have been used previously to determine the shape of slowly oscillating periodic solutions for scalar 
DDEs with constant delay and negative feedback in the singular limit e — )• [27j. They have also been 
instrumental for the analysis of metastable solutions in scalar DDEs with constant delays and monotone 
feedback in the same singular limit [18j . 

For negative feedbacks, the singular limit as e — )• of Eq. ([1]) in the case 7y(0) 7^ can be analyzed 
through the theory developed in |281 [29| [30] that replaces transition layer equations with the so-called "Max- 
Plus" equations. However, given that, as argued in section [3.31 and numerically shown in section [5l such 
systems do not support metastability, we will not dwell any further in this case. Throughout the remainder 
of this section, our focus is on the case 7/(0) = 0, which we henceforth assume to hold. 

For some feedback functions / and state dependent delay function r, solutions of Eq. ([!]) have an 
approximately periodic square-wave shape when e — )• 0, as in the constant delay case r(x,e) = 1, and the 
corresponding "jumps" can be analyzed with the help of transition layer equations. The approximate period 
of metastable oscillations (depending on e) is an essential point in finding transition layer equations. For 
constant delay, at first order, this period is 2 -|- pe for negative feedback ( [27j theorem 3.2), and 1 + ep for 
positive feedback |18j. In section [5] similar asymptotics are shown to hold in the state-dependent delay case 
as well when r]'{0) = and < 77'(0) < -|-oo . In this section we write appropriate transition layer equations 
under various hypotheses on rj, show that their solutions exist, and that these equations can be used to 
characterize metastability, as confirmed by the numerical investigation presented in section [5l 

The case 7/(0) = and 77' (0) = -|-cxd has also been investigated, using r{x, e) = l + e°R{x) with a = ^. We 
found numerically that oscillations are square- wave-like when e is small, their period is either 2-|-e"/9 (negative 
feedback) or 1 -|- e"p (positive feedback), and rescaling time as (/)(t) = x(^), one observes convergence to a 
transition layer profile (see figure [7] in section [5]). However the scaling argument used in the case r/'(0) < -|-oo 
does not apply here, and these transition layer profiles are not analyzed in this section. 

The remainder of this section is organized as follows. First, in section 14. 1[ we show that when < 
77'(0) < -|-oo, appropriately defined transition layer equations can be used to find a symmetry condition that 
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characterizes precisely the cases of metastable oscihatory transients. Details about the existence of transition 
layer solutions, their numerical construction, and illustrating figures can be found in the appendix section 
[B1 The case of ry'(O) = is examined in section [321 One can still write a transition layer problem, but it 
does not depend on the delay function R anymore, leading to incorrect results. To overcome this, we have 
introduced a one-parameter family of auxiliary transition layer problems, and thanks to the analysis of the 
corresponding one-parameter family of transition layer solutions we are able to characterize the cases where 
metastability can occur. Finally, in section 14.31 we discuss a new phenomenon: the possibility of a state 
dependent delay giving rise to metastability in equations that do not exhibit such transients when the delay 
is constant. 



4.1 Case ri{0) = and < r]'{0) < +cx) 

4.1.1 Positive feedback 

Let / be of positive feedback type, and consider Eq. ([1]) where r(x, e) = 1 + r]{e)R{x) with r/(0) = and 
< r/'(0) < +00 (without loss of generality we assume that r]'{0) = 1 in the remainder of this section). 
For such delays, our numerical investigations show that metastable oscillations are approximately 1 + ep 
periodic. 

As for DDEs with constant delays (see section[2]), the jumps of these square- wave-like metastable solutions 
connecting respectively b to —a and vice versa are described by: 



</>±(t) 



= -</.±(t) + /(,/.±(t-i?(0±(t)) + /5±)) 



(6) 



where p~ > and p+ > are unknown constants (drift speeds) and the functions (j)^ satisfy the boundary 
conditions linit-^^oo <P^ {t) = b, lim.t-^oo <P~ {t) = —a , limi_!._oo </>'*' (i) = —a, and liuit^oo (j^^ (t) = b. 

If a solution to equation ([6]) exists, it is called a transition layer solution. In contrast to the constant 
delay-case, this transition layer equation ([6]) is a state-dependent equation, and from a theoretical point of 
view, depending on the values of p and R{x), it may be both advanced and delayed. The numerical method 
used for solving Eq. ([6]) is presented in the Appendix IB. II We have defined an operator T whose (stable 
locally attractive) fixed points are solutions of ([6]). 

In general, the constants p~ and p'^ associated to the decreasing and increasing transition layer solutions 
are different. In the Table [1] we display p~ and p^ solutions of equation ([6|) for r/(e) = e and various choices 
of R{x), for both symmetric and non-symmetric positive feedback / (numerical method details given in the 
Appendix IB. ip . We found that in the state dependent case of equation ([6]), as in the constant delay case, 
oscillatory transients are metastable only when p+ = P- Table [1] also shows that /?"'" = p~ is obtained only 
when the positive feedback function is symmetric and the delay R{x) is even. 



(a) 


R{x) = 


R{x) = X 


R{x) = cos(x) 


R{x) = ix(l + x) 


P^ 


0.824 


0.554 


1.752 


0.732 


P 


0.824 


1.158 


1.752 


1.158 




(b) 


R{x) = 


R{x) = X 


R{x) = cos(x) 


R{x) = ^x{l + x) 


P+ 


1.024 


0.690 


1.916 


0.932 


P 


0.664 


0.994 


1.612 


0.884 



Table 1: Drift speeds p^,p^, solutions of equation ([6]) for positive feedbacks /, 77(e) — e and various choices of delay 
R{x). Metastability occurs only when /7+ ^ p~ ■ Table (a) Symmetric Positive Feedback f{x) = i arctan(5x). Table 
(b) Non-symmetric Positive Feedback f{x) = ^ arctan(5(a; — 0.05)) + ^ arctan(0.25). (see appendix IB.ll for details on 
the numerical method and parameters value used.) 



4.1.2 Negative feedback 



We now discuss the transition layer equation for state dependent delayed negative feedback. Numerically, 
transient oscillations are square-wave-like and they have an approximate period T ~ 2 + eC (see section 
[5]). Likewise the case of DDEs with constant delay, the transition layer equations for the increasing and 
decreasing transition layer solutions of DDEs with state dependent delays are coupled : 









Ri 
R( 



-it)) + p-]) 



(7) 



where (j)'^ is increasing and 4>~ is decreasing on IR, with limt^_oo 4'^{t) = ~0) hmt-i-oo 4'^{t) = ^; hmt^_oo 4'^{t) 
b and limj_!.oo 0~(i) = —a, with </>"*" (0) = (/'"(O) = and p^ are unknown real constants (drift speeds). See 
appendix section |B] for numerical solutions of Eq. ([7]) . 

For negative feedbacks, the symmetry condition supporting metastability is always satisfied when the 
delay is constant [18J, and we show that the same holds when the delay is state dependent. In Table [2] 
we display the drift speeds p~^, p~ that are solutions of Eq. ([7]) for rj{e) = e and various choices of R{x). 
Comparison of Tables [1] and [2] shows that when the feedback function / is symmetric, and R{x) is even, 
p"*" = p~ for both positive and negative feedbacks. This happens because when / is symmetric, and R{x) 
is even, the increasing solutions of the transition layer equations for both positive and negative feedback 
coincide (the same happens for the decreasing solutions). 



(a) 


R{x) = 


R{x) = X 


R{x) = cos(x) 


R{x) = ix(l + x) 


P+ 


0.824 


1.172 


1.744 


1.122 


P 


0.824 


0.434 


1.744 


0.714 




(b) 


R{x) = 


R{x) = X 


R{x) = cos(2;) 


R{x) = ix(l + x) 


P^ 


-0.702 


-0.616 


0.280 


-0.646 


P 


2.574 


2.362 


3.570 


2.476 



Table 2: Drift speeds p'^ , p , solutions of equation ([7]) for negative feedbacks /, 77(e) = e and various choices of delay 



R{x). Metastability occurs regardless of the equality /9+ = p~ . 
(b) Non-symmetric negative feedback f{x) — —^ arctan(5(a; 
on the numerical method and parameters value used.) 



(a) Symmetric negative feedback f{x) = — ^ arctan(5a 
f 0.05)) + iarctan(0.25). (see Appendix [B^ for details 



4.2 Case 7^(0) = and 77'(0) = 

4.2.1 Positive Feedback 

If 7/(0) = and r/'(0) = 0, solutions are approximately 1+ep periodic, and the same time rescaling (p{t) = x{et) 
implies <j){t) = -4>{t) + /(0(t - ^R{(t){t)) + p)) , and as e ^ 



0(t) = -0(t) + /(0(t + p)) 



(8) 



where the decreasing and increasing solutions </)^ must satisfy the boundary conditions lim c, 

t—^—oo 

±(0) 



-(t) 



lim 

t— J- + CXD 



(t) = b, lim (f) {t) = b and lim </> (t) 



-a, with 



t-s> 



t-++oo 



and p = p^ > is an unknown 



constant. This equation ^ is the same one found in the constant delay case, for which the existence of 
decreasing and increasing transition layer solutions has been proven in [18j. In particular, when 7/(0) = 
and 77' (0) = 0, the drift speeds p^ are equal to those of the corresponding constant delay case {R{x) = 0). 
As a consequence, we obtain that for positive feedback /, if 77(0) = and 7y'(0) = 0, metastability cannot 
occur if p+ 7^ p~- However, when p+ = p~, metastability may or may not occur. 



To obtain the symmetry requirement for metastability in this case 77(0) = and t?'(0) = 0, we introduce 
the following 1-parameter family of transition layer equations 



0±(t) = -4>^{t) + f{<t>^[t +pf- xRi^^m), 



(9) 



where, A G M is a real parameter. As previously, p^ are some unknown real constants, and (j) are the 
transition layer solutions, that are expected to depend on A. The case A = reproduces the drift speeds 
p^ of the transition layer equation ([8]) . In figure [1] we display the constants p^ as function of A for the 
symmetric positive feedback function /(x) = | arctan(5x) and various choices of R{x). We found that if 
pj^ = p^ holds only for A = (panels (b) and (d) in figure [1]), metastable DITOs are not observed for 
positive values of e. Oscillatory transients are metastable only when p~^ = p~^ on some non-trivial interval 
A G [0,(5] with 5 > (panels (a) and (c) in figure [T]). So this is the new sufficient condition for the existence 
of metastable oscillatory transients when the positive feedback / is symmetric and 77' (0) = r/(0) = 0. 
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Figure 1: Drift speeds pj, solutions of equation (|9]), for A e [0,0.5], using the symmetric positive feedback function 
f{x) = ^ arctan(5a;) and various choices of R{x): (a) R{x) — (constant delay); (b) R{x) = x\ (c) R{x) = cos(x); (d) 
R{x) — 2^{^ + ^)- See appendix IB. II for details on numerical methods and parameters used. 



4.2.2 Negative Feedback 

For 7/(0) = and ^/'(O) g]0; +oo[, we have seen that the symmetry condition supporting metastability always 
holds when the feedback is negative. The same is true if r/'(0) = 77(0) = 0: regardless of the symmetry of the 
feedback / and of the function R{x) in r(x,e) = 1 + ri{e)R{x), rapidly oscillating transients are metastable 
(see section [5|) . 

4.3 Metastability induced by state dependent delay 

In this section we present a new phenomenon: given a constant DDE that does not display metastability, 
it is possible to add a state dependence to the delay so that the resulting state dependent DDE will exhibit 
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metastability. To this end, we consider equation ([T]) with a non-symmetric positive feedback function /, 
and delay function rx{x, e) = 1 + eXR{x), so that A = corresponds to constant delay DDE which does not 
exhibit metastability. 

We did a numerical investigation using the non-symmetric positive feedback function /(x) = 2 arctan(5(x-f 
0.05)) — ^ arctan(0.25). We have used the following functions R{x): R{x) = x, R{x) = cos(x), R{x) = 
2a;(l -|- x). Solving the transition layer equation ([9]) for each function R we have numerically computed 
the A -families of constants p^ and p~^, the parameter A being varied within the interval [—1.0,1.0] (it 
should be remarked that for large A values the numerical solution of the transition layer equation ^ is 
problematic) . Results are displayed in the figure [2j Metastability will occur for those values of A such that 
p^ = Px- Figure [2]^c) {R{x) = cos(a;)) shows that no solution was found such that p^ — Px^ likewise the 
constant delay case (figure [2)||a)), indicating that introducing a state dependent delay may not make up 
for the lack of symmetry of the feedback function /. Nevertheless, figure [2^b) {R{x) = x) and figure [2]^d) 
{R{x) = ^x{l + x)) show that, in these cases, adding state dependence to the delay has resulted in metasta- 
bility. The solution such that p'^ = p^ is Ac ~ 0.5 in the case R{x) = x (see figure [2]^b)), and Ac ~ 1.1 in 
the case R{x) = ^x{l + x) (see figure [2|^d)). 
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Figure 2: Drift speeds p as function of A for non-symmetric positive feedback function f{x) — ^ arctan(5(a; + 0.05)) — 
^arctan(0.25), and delay function R{x): (a) constant delay case R{x) = 0; (b) R{x) = x; (c) R{x) = cos(a;); and (d) 
R{x) = hx{l + x). (See appendix section IbD for details on numerical methods and parameters value used.) 



5 Numerical simulations of equation ([T]) 

To corroborate the characterization of metastable state-dependent DITOs obtained in sections [3] and HI 
we have carried a numerical investigation of equation ([1]), thus completing the analysis of the transient 
dynamics. 

We shall present the results of numerical solutions of Eq. ([1]), for negative and positive feedback functions 
/, and delay functions of the form r(x,e) = 1 + 7y(e)i?(x). We have used ry(e) = e" with a = (7?(0) 7^ 0), 
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a = 1 (77(0) = 0,?7'(0) > 0) , a = 1 (77(0) = 0,V(0) = 00) and a > 1 (77(0) = 0,V(0) = 0). As for R{x), 
we have considered the following cases: R{x) = x, R{x) = x^, R{x) = |x(l + x), R{x) = sin(x), and 
R{x) = cos(x). 

For positive feedback, with constant delay, metastability requires that the feedback function is an odd 
function of x, so in this case we have used f{x) = ^ arctan(5x). For negative feedback case we have used 
both symmetric f{x) = — ^ arctan(5x) and non-symmetric f{x) = — ^ arctan(5(j; + 0.05)) + ^ arctan(0.25). 
The results are qualitatively the same for these two functions, so we shall only show the results for the 
symmetric negative feedback function. 

In the case 77(0) = and r/'(0) = (a > 1), the observations are qualitatively the same. As a increases, 
the results are closer and closer to those observed in the constant delay case. 

The numerical results were checked using first and second order numerical schemes, using time steps dt = 
5.10^^ and dt = 2.10"^, and with linear interpolation for the state dependent delay function. Simulations 
of solutions of equation ([1]) have also been checked using the RADAR-V package in Fortran. The range 
of e values we investigated is e G [0.01,0.1]. We have used the same initial condition on t G [—2,0] for all 
simulations. Metastability is checked by tracking the zeroes of the solutions. In the positive feedback case, 
we say that transient oscillations end when the last pair of zeros of the solution disappears. In the negative 
feedback case, we say that transient oscillations end when the solution has at most one pair of zeroes in any 
interval of length two (called "slow oscillations" ) . 

In the following, in order to show whether the transient oscillations time T^ is of order exp(-), we plot 
e Vs elog(Te). If y{e) = elog{T^) satisfies y{0) > 0, then T^ = e^^^"'"''^"'^^'', meaning that oscillatory transients 
are metastable. On the other hand, if y(e) = elog{T^) satisfies y(0) = 0, then T^ = e°^~', meaning that 
oscillatory transients are not metastable. 

5.1 Case 77(0) ^ 

In this case, even when e is very small, the oscillations do not last for an exponentially long time. As 
discussed in section 13.31 the lack of metastability in our numerical simulations follows from some of the 
results on a particular class of state-dependent DDEs due to Mallet-Paret and Nussbaum [281 ISHl [30] . As e 
tends to zero, numerical convergence to the limit profile shape described by Mallet-Paret and Nussbaum is 
very clear. 

The results are displayed in the figure El positive feedback on the left panels, negative feedback on the 
right panels. The top panels in the figure [3] exhibit the solution profile for e = 0.1,0.01,0.001, and the 
bottom panels display e Vs elog{T^). 

The solution profile displayed in figure [3]) (top panels) show that when 77(0) 7^ 0, for both positive (top- 
left panel) and negative (top-right panel) feedbacks, the oscillations are not square-wave-like even for very 
small e, in contrast to what happens for constant delay. 

The bottom panels of figure [3] show that when 77(0) ^ 0, oscillatory transient duration grows slowly 
when e converges to zero. For positive feedback (figure [3l^c)), when R{x) = |x(l + x) (curve x), DITOs' 
duration never exceeds a few units of time for e > 0.001; for R{x) = cos(x) (curve +) DITOs' duration 
tends to +00 but it does not grow as e~ when e — ;■ (the function y{e) = elog(T) satisfies y(e) — )■ when 
e -^ 0), and they are not metastable in this sense. For negative feedback (figure [3)^d)) the DITOs' duration 

c 

does not grow as e^ when e — >■ 0, meaning that DITOs are not metastable. Moreover, we can see that in 
the negative feedback case the DITOs' duration depends very little on R{x), the curves e Vs elog(re) being 
almost identical for R{x) = x and R{x) = ^x{l + x) (see figure EJ^d) ) . 

5.2 Case r/(0) = with < r]'{0) < +00 

In this case, in addition to the Hopf bifurcation theorem, Cooke's rescaling argument applies to Eq. ([1]) and 
one expects that the rapidly oscillating solutions have large amplitude when e tends to zero. We numerically 
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(a) f(x) = 0.5*atan(5*x), r(x)=1+cos(x) 
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(d) f(x) = -0.5* atan(5*x) 
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Figure 3: Profile of solutions of equation ([T]) (top panels), and oscillatory transient duration (bottom panels), when 
77(0) ^ 0, with delay r(x, e) = 1 + i?(a;), for positive feedback /(a;) = \ arctan(5a;) (panels on the left) and negative 
feedback j{x) = — iarctan(5a;) (panels on the right). Top panels: solution profile for e = 0.1,0.01,0.001, (a) positive 
feedback with R(x) — cos(x), (b) negative feedback with R{x) — x. (c) e Vs e\og{Te) for positive feedback, and delays 
R{x) = cos(a;) (curve +), R{x) — ^x{l + x) (curve x); (d) e Vs elog(Tc) for negative feedback, and delays R{x) — x 
(curve +), R{x) — ^x{l + x) (curve x). 

observed that metastable oscillations are almost 1 + ep periodic (the zeroes drift-speed p depends on the 
feedback function / and the delay function r), and that the constants p are coherent with the corresponding 
constants in section 14.11 This period estimate is crucial for obtaining the transition layer equation for the 
state dependent DDE (see section|H equation ([6|) ). Here we have taken ry(e) = e, so that r{x, e) = l + eR{x). 
Metastability is expected for any delay function R in the negative feedback case (symmetric and non- 
symmetric), while for the symmetric positive feedback case the delay function R must be even so as to 
guarantee the symmetry condition. 

Figure m displays the numerical results for Eq. ([1]) when 77(e) = e, for positive feedback (panels on the 
left) and negative feedback (panels on the right). The top and middle panels of figure H] display the DITOs' 
profiles for R{x) = cos(x) and a few e values. The bottom panels of figure [4] display the transient duration (e 
Vs elog(Te)) for R{x) = (constant delay), R{x) = cos(x), R{x) = x, R{x) =, R{x) = |x(l + x), R{x) = x. 

The top panels of figure H] show that for both, positive and negative feedback cases, the oscillatory 
solutions have a square-wave-like shape when e goes to zero, likewise the constant delay case. Figure H^a) 
for positive feedback shows that, as e ^- 0, the square- wave-like solution has period 1 + ce (at first order). 
Figure mb) for negative feedback shows that, as e — )• 0, the square-wave-like solution has period 2 + c'e (at 
first order). 

The middle panels in figure S] display the square- wave- like oscillation after rescaling the time, s = -. The 
time-rescaled profiles displayed in the middle panels of figure U] show convergence to a limit profile when e 
converges to zero, for both positive and negative feedback, in agreement with the results in section [3l As 
explained in section [3l after the appropriate time rescaling s = -, the square wave shape converges to some 
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(a) f(x) = 0.5*atan(5*x), r(x,e)=1+ ecos(x) 



(b) f(x) = -0.5*atan(5*x), r(x,E)=1+ ecos(x) 
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Figure 4: Profile of solutions of ([TJ (top and middle panels), and oscillatory transient duration (bottom panels), when 
77(e) — e, for positive feedback f{x) — \ arctan(5a;) (panels on the left), and negative feedback f{x) = — i arctan(5a:) 
(panels on the right). Top panels: solution profile for R{x) = cos(x) and decreasing e. Middle panels: descending 
and ascending jumps obtained by the time scale s = - of the solutions profile for R{x) — cos(a;), varying e: (c) 
positive feedback; (d) negative feedback. The curves for e — 0.03 and e = 0.001 are too close to be distinguished. 
Bottom panels: (e) transient duration for positive feedback with delays R{x) = 0, R{x) = cos(a;), R{x) = x, 
R{x) = \x{\ + a;), R{x) — x; (f) transient duration for negative feedback with delays R{x) — 0, R{x) = cos(.t), 
R{x) = X, R{x) — ^x{l + x). 

limit profile, which is solution of the transition layer problem ([6]). Convergence to transition layer profiles 
occurs for both increasing and decreasing rescaled jumps, for both positive and negative feedbacks. 

An approximation of the unknown constant p = p in problem ([6]) can be obtained from the numerical 
approximate period of metastable oscillations : T ~ 1 + ep. We have checked that the limit profiles agree 
with the solutions of the transition layer equations obtained as attractive fixed point of an appropriate 
operator T, up to numerical error of order 0{dt/e), where dt is the discretization parameter (see appendix 
section [B]). Due to regularity properties of the operator T, this implies that the corresponding constants p 
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also agree at the same order. 

The figure Sf shows that in the negative feedback case, metastable oscihatory patterns are observed 
regardless of the delay function R. For negative feedback with constant delay {R{x) = 0), the curve + in 
figurelD^f) shows that T(e) ^ exp(-), for some non-zero constant c, as expected. This same panel (f) shows 
that, for state dependent delay, elog(T(e)) r^ ci + C2e + o(e), implying that T(e) = exp(^) exp(c2)(l + o(l)). 
The value of constants ci and C2 depends on the delay function R, but in all cases ci > implying that the 
state-dependent DITOs' are metastable. 

State-dependent DITOs are metastable if and only if the function y[e) = elog(T) has a non-zero limit 
when e — >■ 0. Figure |4^ shows that in the positive feedback case, the DITOs' are metastable only for 
odd functions / and even delay functions R{x). The curve + in the figure IH^e) shows that, for constant 
delay, T(e) '^ exp(|), for some non-zero constant c, as expected. In the case of state dependent delay, 
when the delay function is even, the curves for R{x) = cos(x), R{x) = x^ in the figure Sl^e) show that 
elog(r(e)) Rr: ci + C2e + o(e), meaning that T = exp(^) exp(c2)(H-o(l)) so that the state-dependent DITOs' 
are metastable. In contrast, when R{x) is not even (cases R{x) = sin(3;) , R{x) = x + ^x^), figure Hl^e) 
shows that elog(r(e)) — > so that r(e) is not of order exp(-), implying that the state-dependent DITOs' 

are not metastable. 

5.3 Metastability induced by state dependent delay 

As already emphasized, for positive feedback /, DDE ([1]) with constant delay exhibits metastability only 
if / is odd. The constants p'^ and p~ that solve the TLE problem are equal in the case of odd positive 
feedback / with constant delay, implying metastability. When the positive feedback / is not symmetric, 
the condition p"*" = p~ does not hold for constant delays, and DITOs are not metastable. Nevertheless, 
in subsection 14.31 we have shown that for non-symmetric positive feedback /, and state dependent delay 
r(x,e) = 1 + XeR{x), for a critical value A = Ac, the transition layer equations ([9]) do have solutions such 
p'^ = p, provided R{x) is not even (see figure [2]). Therefore, metastability has been induced by introducing 
the appropriate state-dependence to DDE ([1]) with non-symmetric positive feedback. 

In this section, we present the numerical solutions of equation ([T]) for the same case analyzed in section [43l 
f{x) = ^arctan(5(x + 0.05)) — ^ arctan(0.25), with state-dependent delay r{x,e) = 1 + XeR{x), R{x) not 
even. Panels (a) and (c) of figure [5] display the results for the case R{x) = \ex{l + x) (corresponding to 
figure [2jd)); panels (b) and (d) display the results for R{x) = x (corresponding to figure [2jb)). 

The fronts and the values of p^ and p^ change continuously with A. In other words, simulations over 
fixed durations of DDE solutions for A in the vicinity of Ac are similar because for A close to Ac, the diff^erence 
/9+ — p~ is small and the DITOs are long lasting. This similarity not withstanding, the transient regime 
durations scale differently with e at Ac and nearby values A. Only at Ac the system becomes metastable in 
the sense that the DITOs last for exponentially long times. This difference is illustrated in the figures that 
show transient regime duration at fixed A for various e and the reverse, i.e. at fixed e for various values of A. 

The panels (a) and (b) of figure [5] display e Vs l/log(T(e), for three values of A, and we can see that 
when the parameter A is larger or smaller than the critical value Ac, the curves y{e) = l/log(T£) are very 
steep when e ^ 0, indicating that DITOs are not metastable in those cases. In contrast, when A = Ac, 
the curves y(e) have a bounded slope as e ^' 0, indicating that DITOs are metastable. We have obtained 
Ac ~ 1.05 for the case R{x) = ^ex(l + x) (figure [5)^a)), and Ac ~ 0.55 for the case R{x) = x (figure [5jb)). 
The difference between these values of Ac, and those values of Ac found in subsection 14. 3^ is smaller than 
10~^, which is of the order of numerical precision for these parameters. 

The panels (c) and (d) of figure [5] display A Vs elog(Te), for three values of e, and we can see that when 
e is fixed, there is a unique value A*^ such that DITOs' duration is maximal, and it converges to the critical 
value Ac as e decreases. 

The results in this subsection agree with the results in subsection 14. 3( confirming the possibility of state- 
dependence of the delay inducing metastability when the constant delay case does not exhibit metastability. 
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(a) r(x,e) = 1 + 0.5 A,sx(1+x) 








(b)r(x,E) = 1 +XE-X. 
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Figure 5: DITO duration for positive non-symmetric feedback f{x) = i arctan(5(a:; + 0.05)) — ^ arctan(0.25), with 
state dependent delay r{x, e) = 1 + XeR{x). Panels (a) and (c): R{x) — ^x{l + x); panels (b) and (d): R{x) — x. 

5.4 Case 77(0) = with ^(0) = 

When 7?'(0) = 0, likewise the case < r/'(0) < 00, oscillations have a square wave shape, a "period" 
T = 1 + ep, and a rescaled limit transition-layer profile (figures not shown). The corresponding transition 
layer equation ([8]) is independent of the function r/, and it is the same transition layer equation as for the 
constant-delay case r{x,e) = 1. 

For given feedback / and delay R{x), if metastability occurs when < 77' (0) < -|-cx3, our numerical 
investigation indicates that it will also occur for r/'(0) = 0. We have computed the transient duration for 
the symmetric positive feedback f{x) = \ arctan(5x), with delay r(x, e) = 1 -|- r]{e)R{x), r/(e) = e", q > 1. 

Figure [6] (a) displays e Vs elog(Te) for the state dependent delay R{x) = |x(l -|- x), and a = 1.5,2,3. 
It shows that elog(T) converges to zero when e — > 0, regardless of a, implying that for symmetric positive 
feedback the DITOs are not metastable if R{x) is not even. 

Figured] (b) displays e Vs elog(r£) for the state dependent delay R{x) = x^ (even function), and a = 2,3. 
It shows that elog(T) does not converge to zero when e — ?> 0, implying that for symmetric positive feedback 
the DITOs are metastable if R{x) is even. From Figured] (b) we can also see that as a increases, the duration 
T of transient oscillations is closer and closer to the duration of the transient for the constant delay case. 

For negative feedback the same result holds: metastable DITOs' for the case < r/'(0) < -l-cxo implies 
metastable DITOs' for the case 7?'(0) = . 

5.5 Case r/(0) = with /^'(O) = +00 

When r/(0) = 0, with rj(^) = -|-oo oscillations have a square-wave-like shape and an approximate period 
r = 1 + e"/> when e — ?> (figures not shown). The usual time rescaling s = - does not give converging 
profiles, but the scaling s = -^ does. As illustration, we display in figure [7] the results for the negative 
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(a) f(x) = 0.5*atan(5*x), Non-even delays 
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Figure 6: Oscillatory transient duration for symmetric positive feedback f{x) = ^arctan(5x), with constant delay 
(curve +), and state dependent delay r{x,e) = 1 + r/(e)i?(a;), 77(e) = e", a > 1: (a) R{x) = ^x{l +x), a = 1.5,2,3; 



(b) R{x) =a;2, a = 2,3 



feedback f{x) = —^ arctan(5x), with r = 1 + e^'^ 003(2;). At the time scale te ^'^, we observe convergence 
of oscillating solutions jumps to a limit transition layer profile when e — > 0. 

f(x)= -0.5*atan(5*x), r(x,E) = 1+ e"^cos(x) 




Figure 7: Convergence of rescaled jumps of oscillating solutions in the limit e -> 0, when ry(0) — with ri'{0) = +cx3. 
We used negative feedback f{x) = — i arctan(5a;), with r = 1 + e^^^ cos(x). 

Here neither Mallet-Paret and Nussbaum's nor Cooke's argument apply, so that we have no indication on 
whether rapidly oscillating periodic solutions (which are expected because of Eichmann's Hopf bifurcation 
theorem) have large or sinall amplitudes. Numerical simulations have shown that oscillatory transients can 
last for a very long time, but due to numerical difficulties it is not clear whether this transient duration is 
exponential, hence we cannot say whether these long lasting oscillatory transients are indeed metastable or 
not. 

6 Discussion and conclusion 

At present time, no global geometric characterization of the organization of the phase portrait of scalar 
state dependent DDEs is at hand. Even basic results, such as the Hopf bifurcation theorem, have only 
recently been established. However, based on current knowledge and systematic numerical explorations, it 
is possible as done in section [3] to conjecture that the phase portraits of DDEs with state dependent delays 
have the same geometrical organization as those of constant delays. Furthermore, it is possible to discuss the 
putative occurrence of metastability based upon available information regarding rapidly oscillating periodic 
solutions, and educated guesses supported by numerical investigations. 

In this paper we have shown that metastable oscillating solutions can exist in singularly perturbed DDEs 
with state dependent delays of type ([I]). Based both on mathematical and numerical results, we have been 
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able to link the properties of Eq. ([T]), and its solutions, to the occurrence of oscillatory metastable transients. 

Such metastable transients have never been observed when the delay r = r(x) does not depend on the 
singular parameter e. The main reason for the short duration of such DITOs' is the decreasing size of 
rapidly oscillating periodic oscillations. Considering delays of the form r(a;,e) = 1 + r]{e)R{x), we found 
that the scaling parameter value r/(e) is crucial for the existence of metastable transients. When rj{0) ^ 
exponentially long lasting oscillatory transients were never observed, even for very small e. When r](0) = 
and //'(O) = oo long lasting oscillatory transients can be observed, but due to numerical difficulties one 
cannot conclude whether these transients are indeed metastable or not. For 77(e) = e", the larger a is, the 
longer transient oscillations will last. When r](0) = and < ^^'(O) < cx) metastable oscillatory transients 
can always be observed for negative feedback, while for positive feedback some symmetry condition must 
be satisfied by the feedback function / and the delay function R. 

There are two main tools to analyze metastability phenomena. With a geometric approach, one can look 
for the existence of a global attractor containing a "cascade" of unstable periodic orbits and heteroclinic 
connections between them; and with a dynamical approach one can investigate transition layer equations 
that describe the asymptotic shape of the oscillations when e converges to zero. 

We have shown that, with few hypotheses on the functions / and r, a Hopf bifurcation theorem of 
Eichmann applies to Eq. ([T|). This implies the existence of a sequence of Hopf bifurcations as e converges 
to zero, meaning that a global attractor with the previously described structure might exist in many cases. 
Nevertheless, metastable oscillations cannot be observed in general. Our numerical investigation has shown 
that when the delay function r does not converge to a constant as e converges to zero (see figure [3] in 
subsection 15. ip . the amplitude of periodic solutions has to converge to zero as the period converge to zero, 
and metastable oscillations were not observed. This suggests that not only rapidly oscillating but also large 
amplitude periodic solutions are needed to support metastability. 

Furthermore, even when the delay function r converges to a constant as e tends to zero (section \5\ cases 
77(0) = 0), the state dependent DITOs' are metastable only in those cases where there exist transition layer 
equations similar to those found for DDEs of type ([T|) with constant delay (a > 1, subsections 15. 4| and l5.2p . 
This suggests that metastable state dependent DITOs' cannot exist unless the oscillations have a limiting 
shape that are solution to a transition layer problem of the form ([6|), as e converges to zero. 

Our numerical analysis of the transient oscillations in DDE ([T]) induced by state dependent delay of the 
form r{x) = l+r]{e)R{x), has revealed that for negative feedback /, metastable state dependent DITOs' exist 
in the same way as for the constant delay case, while for positive feedback /, metastable state dependent 
DITOs' exist if and only if / and R satisfy some symmetry conditions. If the positive feedback / is an 
odd function (symmetry requirement for metastability in the case of constant delay), then metastable state 
dependent DITOs' exist only if the delay R{x) is an even function. We also have showed that by adding 
state dependence to the delay, it is possible to obtain metastability for positive feedback / for which the 
constant delay transient oscillations are not metastable (see subsections 14.31 and 15. 3p . 

An important contribution of this work has been the introduction of a novel class of transition layer 
equations associated with state dependent delays. Prom our numerical investigations, we claim that these 
equations capture two essential aspects of the dynamics of DDEs with state dependent delays. The first is 
the shape of the oscillations as the parameter e becomes small. The second is the drift of the oscillations. 
Our focus has been on monotone feedbacks, nevertheless the transition layer equations we have introduced 
remain valid for non monotone feedbacks as well. Our paper paves the way to investigate the dynamics of 
DDEs with such feedbacks and state dependent delays, through the novel transition layer equations. 

A Proof of theorem [1] 

The proof of theorem [1] is based on a "Hopf bifurcation theorem" proven by Eichmann (Eichmann's PhD the- 
sis [12j p81). PoreG]0;l[andy G Ci([-M, 0],E), let 5 be the function ^(e,?/) = i (-y(O) + /(y(-r(e,y)))). 
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such that Eq. ([T]) is equivalent to 

x'{t)=g{e,xt), (10) 

with the notation Xj(s) = x{t + s) for all s € [—M, 0]. 

The Hopf bifurcation theorem of Eichniann has three first order derivatives hypotheses {Hi) {H2) and 
(i^a), three second order derivatives hypotheses (^^4) {H^) (Hq) and three spectral hypothesis (Li), (L2) 
and (I/s). For convenience of the reader we state these hypotheses here before proving that they are satisfied 
by equations ([I]) and ([2]). 

Suppose that there is an open subset U of C^{[—M,0],M.) such that 

• (^1) the mapping 51 : ]0, l[xC/ ^^ M is continuously differentiable, 

• (-^2) for any (e,x) g]0, l[xC/ the second partial derivative Dxg{€,x) can be extended to a linear 
continuous map 

D,g{e,x): C'>{[-M,0],R) ^R ; 

• (^3) the (extended) mapping 

]0,l[x?7x C°([-M,0],M) ^M 

(e,2;,y) ^ D^g{e,x)y, 

is continuous; 

• (^4) the mapping g : ]0, l[x ([/ n C^) — > M is twice continuously differentiable; 

• (-ffs) for any (e,x) g]0, 1[xU n C^ the second order partial derivative D'^g{e,x) can be extended to a 
bilinear continuous map 



Dlg{e,x) : C\[-M,0],R) x C\[-M, 0],R) -^ R ; 

• {Hq) the (extended) mapping 

]0, l[xU xC\[-M, 0],R) X C\[-M, 0],R) ^R 

{e,x,y,z) ^ Dlg{e,x){y,z) 

is continuous, and 

]0,l[xU xC^{[-M,0],R) ^£{C^,R) 

(e,x,y) ^ Dlg{e,x){y,-), 

is continuous (where C{C'^,R) is the space of linear functionals from C^([— Af, 0],M) to R). 

Suppose that g{e,0) = 0, for any e €]0, 1[, and let ^(e) be the generator of the strongly continuous 
semigroup on C^{[—AI,0],R) generated by the linearized equation 

y't = D,gie,0)yt. 

Suppose that there is a e* g]0, 1[ and some open interval ]e* — ??, e* + ??[c]0, 1[ such that 

• (Li) for any e G.]e* — r], e* + r][, there is a simple eigenvalue A(e) of A(e), such that the mapping e 1-^ A(e) 

isC\]e*-v,e* + 7i[,C), 



(L2) the eigenvalue A(e) crosses the imaginary axis at e* : K(A(e*)) = and 3(A(e*)) = (Jq > 0, and 

(L3) and for any A; € Z — { — 1, 1}, i' = ikujQ is not an eigenvalue of A{e*). 
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Now we prove that theorem [T] is a consequence of Eichmann's theorem. To this end we start by proving 
the following Propositions [3] and [H 



Proposition 3. Let U be an open subset o/ C^([-M, 0],R). Suppose f : R ^ R is C^ and r :]0; l[xU — )• 

0],M) -^ 

, and that (e, xt, ht) — )• ^{^, xt).ht is continuous ]0; l[xU x 



0,M] is C^. Suppose that for any {e,xt) G]0;1[x?7, |^ :]0; 1[xC^([-M,0],M) -^ R can be extended to a 



continuous linear map ]0; l[xC^([— M, 0],M) - 
CO([-M, 0],M) ^R. 

Then the regularity hypotheses Hi, H2 and H^ of Eichmann's H op f Bifurcation Theorem hold. 

Here, ^ and ^ denote Frechet derivatives. 



Proof. One finds 



do — 1 1 t dr 

■^(e,2/) = —g{e,y) - - {-g^{f^,y)y'i-r{e,y))f'{y{-r{e,y))) 



-—{e,x).h = - 
oy e 



-h{0) + (-x'i-r{e, x))^{e, x).h + /i(-r(e, x))\ f'{x{-r{e, x))) 



This shows that §f :]0, l[x[/ ^ M and || :]0, l[x[/ ^ M are C°, so that g G C\ and Hi holds. Using 
the hypotheses made on r in the theorem [H the formula shows that H2 and H^ also hold. D 

Proposition 4. Suppose that the hypotheses of Proposition holds. Suppose additionally that / : M — )• M 
is C^ and that r :]0,l[x ([/ n C2([-M,0],M)) -^ [0,M] C M is C^. Suppose that for any (e,x) G]0,l[xf/, 

the function jT^{e,x) has a bilinear continuous extension to C^([— M, 0],M) x C"'^([— M, 0],M) — > M, which 
depends continuously on (e, x) €]0, l[x(C/ n C ). ^nd suppose that {e,x,h) \-^ ^{e,x){h,-) is continuous 
]0,l[x(C/nC2) xC^ -^C{C^,R). 

Then regularity hypotheses H4, H^ and Hq of Eichmann's Hopf Bifurcation Theorem are true. 



Proof. This proposition is a consequence of the following formulas for second order derivatives of function 

g- 

d ( dg 



^ .e-+5)(e,x) 



dh 



ae2 



(e, x)x'(-r(e, x))/'(x(-r(e, x))) 



-x"{-r{e,x)){—{e,x)\ /'(x(-r(e,x))) - f x'(-r(e,x)) — (e,x) j /"(x(-r(e,x))) 



and 



-^y{'% + a]{e.x).h 



8h 
dyde' 



{e,x).h 1 x'(-r(e,x))/'(x(-r(e,x))) 



-■^(e,2;)a;"(-r(e,x)) ( — (e,x)./i) /'(x(-r(e,x))) 



5r 



+— (e, x)x'(-r(e, x))/"(x(-r(e, x))) 



-x'(-r(e, x)) f — (e, x)./i ) + h{-r{e, x)) 



and 



9i V 9y ^^' ^) ■ ^ ) = /'(^(-^(«' ^))) 



OT ( OT 

x"(-r(e,x)) — (e,x) ( — (e,x)./i 



+/'(x(-r(e,x))) 



9r 



-x'(-r(e,x))— — -(e,x)./i - /i'(-r(e,x)) — (e,x) 



-x'(— r(e,x))— -(e,x)./i + h(—r{e,x)) 
oy 



/"(x(-r(e, x)))x'(-r(e, 2;))^(e, 2;) 
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and 



- — (e, x).h + h{-r{e, x)) j ( -^(e, x).k + k{-r{e, x)) ] x'{-r{e, x))f"{x{-r{e, x))) 



+/'(x(-r(e,x))) 



x"(-r(e,x)) ( ■^{e,x).kj l^-^{e,x).h 



-f\x{-r{e,x))) 



d T 3t 

-x' {-r{e, x))-—^{e, x).{h, k) - h'{-r{e,x)) — {e,x).k 



If r and / are C^, then one can check that g :]0, l[x ([/ n C2([-M, 0], M)) ^ M is C^ and this is H^. If 

|^(e,rc) has a continuous extension to C^{[-M,(}\,R) x Ci([-M,0],M) ^^ M, then so does |^(e,x) and this 

is H^. Since the extension of ^{e,x) and ^-^(e, x) are continuous in {e,x) the first part of Hq is satisfied. 

Since r :]0, l[xU — )■ M is C , and due to the last hypothesis on ^ in the proposition, the last requirement 
of Hq holds too. D 

The hypotheses of theoremU /isC2(M,M),r :]0, l[xC°([-Af,0],M) ^MisC\andr ■.]0,l[xC^{[-M,0],l 
M is C^, imply that the regularity hypotheses of Proposition [3] and Proposition H] are verified. 

We now turn to the spectral hypotheses of Eichmann's Hopf-bifurcation theorem. We consider the 
equilibrium x* = 0, which satisfies g{€,0) = for all e g]0, 1[, and the linearized equation at x* : y'{t) = 
|(e,0)2/i withyoGCi([-M,0],M), ie 

ey'{t) = -y{t) + /'(O)y(t - tq) 

where rg = r(0) = r(x*) = 1. The characteristic equation is 

1 + ea = /'(0)e-"cos(/3) 
e/3 = -/'(0)e~"sin^, 

This is the same characteristic equation for the constant-delay equation, and a standard argument shows 
that there exists a sequence e^ — y such that for each e = e^ the characteristic equation has a single pair 

fc— >+oo 

of solutions on the imaginary axis A = ±i/3fc, /3fc > 0. This implies in particular that (L3) of |12] is satisfied. 
To check that for each k the eigenvalue A^ can be tracked in a neighborhood of e = e^, we use an implicit 
function theorem. Considering G :]0, l[xM x M^ defined by G(e,a,/3) = {1 + ea — /'(0)e~" cos(/3); e/3 + 
/'(0)e~" sin(/3)), we have G(e,a,/3) = if and only if A = a + i/3 is a characteristic root. We have 
G(efc,0,/3fc) = (for any A;), and the derivative of G with respect to a and j3 at (efc,0, /3fc) is 

/ l + efc + /'(0)cos(/3,,) /'(O) sin(/3fc) \^( 2 + e^ /'(0)sin(/3fc) A 

V -/'(0)sin(/3fc) e, + /'(o)cos(/3fc) y V -/'(0)sin(/3fc) 1 + e^ ^ 



which has a positive determinant and is invertible. Thus (for any k) there is an open interval I^ containing 
efc and a C^ function v^ : I^ ^ C such that for all e G 1^, ffc(e) is a characteristic root and I'ki^k) = ^k- 
Furthermore, computing the e derivative of ffc at e = e^ and using the relations for Afc, one finds 

(1 + efc)a'(efc) = efe/3fc/3fc(efc) 

(l + efc)/3^ = -il + eka'^iek))/3k, 

which gives either /3^(efc) = and a'f,{ek) = 7^, or /3^(efc) 7^ and a;'^(efc) = l^-^-ill^^ so that one always has 
a^(efc) ^ 0. This means both (Li) and (L2) of Eichmann's spectral hypotheses are satisfied, which finishes 
the proof of theorem [TJ 
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B A numerical method for solving transition layer equations. 

In this section we present a method for solving numerically the transition layer equations ([6]), ([7]), ([8]) and 
dl]). This section is divided into two parts: negative and positive feedback. Some details are only provided 
for the positive feedback case, since for the negative feedback case they are similar. 

B.l Positive feedback 

Existence of transition layer solutions is usually proven with the help of a fixed point theorem. Suppose that 
/ is smooth, has three fixed points x = —a, x = and x = b, and that / is increasing on the interval [—a, b] 
(positive feedback, with no symmetry hypothesis). Consider the set C^ of functions (p G C^(M,M) such that 
(/>(0) = 0, lim (/)(x) = —a, lim (p{x) = b, ||(/)'||ioo < sup |/(3;)| = max{a, 6}, and (j) is strictly increasing 

^-^~°° x^+oo xe[-a;b] 

on R. On that space we define the operator T that to G C^ associates the function T(j) = ip given by 

^(t) = e-* / e7(0[s + P- RiHsWds, 

J —oo 

where p G [0; +oo[ is the only constant such that tpi^O) = f_^ e*/((/>[s + p — i?(</>(s))])(is = 0. It can be shown 
that ip G C^. Then ip = T(j) is the unique solution of equation 

m = -^{t) + f{4>[t - Rim) + p]) (11) 



such that lim ip{x) = —a, lim ^(x) = b and V'(O) = 0. So an increasing solution of the transition layer 

x— >— oo X— >+oo 

equation ([6]) (case r/(0) = and r?'(0) = 1), is a fixed point of the operator T in the set C+. In the following 
we show how to numerically solve the fixed point problem T</>''' = (j)^. A solution to the problem T(/>~ = (j)~ 
is obtained in the same way. 

Let 4> he a smooth increasing function such that lim (/)(t) = —a and lim (p{t) = b. Then, for any 

i— >— oo t— >+oo 

7 £ M, the function 

-t 



V^(t) = e-* / e'f{4>[s + 7 - RiHsWds 

J —oo 

is also smooth, increasing, and satisfies lim ^(t) = —a and lim ip{t) = b. Since, for any t € M, the 

t— >— oo t— >+oo 

map 7 i-7> ip-y{t) is increasing, there is a unique 7 = /o such that ippiO) = and so Vp = 'T4'- A numerical 
approximation to the map (j) ^ ipp is the following. Let L > and (j)^ be the function that coincides 
with (p on the interval [— L,L] and such that (p (t) = —a for t < —L and (/) {t) = b for t > L. Clearly 
suPigR I'P^i't) ~ 0(01 can be made arbitrarily small if L is chosen sufficiently large. So, we fix L > and 
choose an initial 0o that satisfies 0o(O — ~^ lo'^ ^ < ~^ ^'^d (poit) = & for t > L. Then for a given 7 we use 
a first order Euler method to solve the equation 

v^;(t) = -^^{t) + f{4>[t+j- Rium) , 

with the initial condition ip^[—L) = —a and time step dt. Knowing that 7 > /? is equivalent to ^^(O) > 0, we 
can use a shooting method to find 71 and 72 such that p — (it < 71 < /) < 72 < /o + dt. Defining (j)i{t) = ^^^ (t) 
for t < L and (j)i{t) = b for x > L we obtain the approximation T4>q ~ 4>i. This procedure can be iterated 
(pn+i = Tcpn, n = 1,2,... hoping that it converges to a fixed point cp'^ . The convergence of this sequence 
was numerically verified, typically iterating (pn+i = Tcpn up to n = 20 and checking that 

|0„(t)-0„_i(t)| <dt. 

A sample of our results are shown in table [1] (a) (section 14. l.ip and in figure [8j 
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In table [T] (a) the constants p^ were computed for f{x) 
The initial condition used to obtain table [1] (a) was (p^ (t) = 



iarctan(5a;), L = 100.0, and dt = 0.001. 



a&(l-e~?') 



a+be 3 






for the increasing transition layer 



solution, and (p^ (t) = — — %j^ for the decreasing transition layer solutions. 

a+feeS 

In figure [5] (a) several iterates (/>„ are shown to converge to a limit profile (f)' for f{x) = i arctan(5x), 
with R{x) = cos{x), and initial condition (t>o{t) = ab ^~^^i . In that case the feedback / is symmetric, so that 
the operator T is also symmetric, and if we consider an increasing initial profile (J)q = —4>o, the corresponding 
sequence is (pn = T'^{—<Po) = —Tcpo- In particular, figure [8] illustrates the convergence and shapes of iterates 
for both increasing and decreasing profiles. In figure [8] (b) the limit profile (j)io{t) = T^^cpoit) ^ 4'~{i) is 
shown for R{x) = 2x, R{x) = — x^, and R{x) = \x{l + x), using the initial condition 4>o{t) = ci^ a+be^ - 
For R{x) = X, the numerical convergence of the sequence (pn+i = Tcpn was successfully tested for the 
following initial functions: <^o(i) = ~if[) ^-^^ 4>o{t) = CLb ^~^^i (1 + 0.4cos(t)). All these results also hold for 
a non-symmetric positive feedback as well. 



(b) 
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Figure 8: Decreasing transition layer solutions of equation ([6|) (?/(0) — with r]'{Q) — 1) with f{x) ~ i arctan(5a;) 
(Positive feedback case). In (a): R{x) — cos{x) and pait) = oh ^~^^t ■ In (c): 0io(i) — T^'^poi't) is displayed (with 
poit) = ah Iv^^t ) , for functions R{x) = 2x, R{x) = —x^, and R{x) = jx{l + x). 

For DDEs with state dependent delays, such as transition layer equations ([6]), the existence of solutions 
(j) to the Cauchy problem is known under classical hypotheses, in particular on the delay function R, that 
ensures that t i— >■ t — R{(j){t)) is not decreasing (see [19] for example). We mention that in several numerical 
examples, when R is too large, the maps t ^^ t — R{4>n{t)) are not monotone, the iterative sequences 
4)n = T'^(pQ does not converge, and transition layer solutions seem not to exist (see section 14. 2p 

In the case r/(0) = ?/'(0) = 0, the transition layer equation ([9]) is associated to the operator 



TxHt) 



-t 



e'f{(P[s + p+-\R{cPismds, 



and the same numerical method presented above can be used to compute its fixed point. For instance, 
figure [Tj (section 14. 2. ip shows the constants p that were computed using various choices of R and: /(x) = 

^ arctan(5a;), L = 100.0, dt = 0.001, n = 20, and the initial condition (pQit) = "" ^ ~^ ^ ' for the increasing 

a+fee~3 

2 J 

— — %—^ for the decreasing transition layer. 

Finally, in figure [2] (in section H3|) . illustrating that metastability may be induced by state dependent 
delay in the case of non-symmetric positive feedback (that does not exhibit metastability for constant 
delay), the values of p^ were computed using L = 100.0, dt = 0.001, n = 20, and the initial conditions 



transition layer solution and 0g (t) 



a+be 3 



a+be3' 
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B.2 Negative feedback 

For a negative feedback /, one can repeat the procedure above and define the operator T that has as fixed 
points the solutions to the transition layer equation ([7]), namely 






where (p~^ (resp. (j) ) ^^^ increasing (resp. decreasing) smooth functions with lim 0^ = —a, lim cp'^ = b 
and lim (j)~ = b, lim (p^ = —a, and p^ are the only constants such that 4>'^{0) = 4>^{0) = 0. Since the 

t—>-—oo t— >+oo 

function / is decreasing, the operator T now maps an increasing function to a decreasing one and vice versa. 
As in the case of positive feedback, if the function (j) is strictly monotone the function 7 1— > ^-f{t) is strictly 
monotone too, so that the constants p above are well defined and can be computed by a shooting method. 
We choose a smooth increasing initial function (J)q with lim (/>"'" = —a, lim (j)^ = b and the subsequences 

t— >— 00 t— >+oo 

(j)2n and p2n (resp. 02n+i and P2n+i) converge to the transition layer solution cp'^ and the constant p"*" 
(resp. (j)~ and p~ ). Numerically, as in the positive feedback case, we use discretization of step dt on an 
interval [—L^L], the functions ^^ are computed using a first order Euler scheme of step dt, the constant p 
are approximated with a precision dt, and the convergence after n iterations of T is checked similarly : 

|0„(t)-0„_2(t)l < dt, 

\(t)i-i{t) - h-^{t)\ < dt. 

In table [2] (section I4.1.2p the constants p^ were computed using L = 100.0, dt = 0.001, n = 40, and 

a+be 3* 
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